NASA Technical Memorandum 102879 


Estimates of Thermochemical 
Relaxation Lengths Behind 
Normal Shock Waves Relevant to 
Manned Lunar and Mars Return 
Missions, the Aeroassist Flight 
Experiment, and Mars Entry 

John T. Howe 


(NASA— TM- 102879) ESTIMATES OF N91-21471 

THERMOCHEMICAL RELAXATION LENGTHS BEHIND 
NORMAL SHOCK WAVES RELEVANT TO MANNED LUNAR 

AND MARS RETURN MISSIONS, THE AEROASSIST Unci as 

FLIGHT EXPERIMENT, AND MARS ENTRY (NASA) G3/34 0008215 


February 1991 


NASA 

National Aeronautics and 
Space Administration 






NASA Technical Memorandum 102879 


Estimates of Thermochemical 
Relaxation Lengths Behind 
Normal Shock Waves Relevant to 
Manned Lunar and Mars Return 
Missions, the Aeroassist Flight 
Experiment, and Mars Entry 

John T. Howe, Ames Research Center, Moffett Field, California 


February 1991 


NASA 

National Aeronautics and 
Space Administration 

Ames Research Center 

Moffett Reid, California 94035-1000 




SUMMARY 


Thermochemical relaxation distances behind the strong normal shock waves associated with 
vehicles that enter the Earth’s atmosphere upon returning from a manned lunar or Mars mission are 
estimated. The relaxation distances for a Mars entry are estimated as well, in order to highlight the 
extent of the relaxation phenomena early in currently envisioned space exploration studies. The 
thermochemical relaxation length for the Aeroassist Flight Experiment is also considered. These 
estimates provide an indication as to whether finite relaxation needs to be considered in subsequent 
detailed analyses. For the Mars entry, relaxation phenomena that are fully coupled to the flow field 
equations are used. The relaxation-distance estimates can be scaled to flight conditions other than 
those discussed. 
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universal gas constant 

number of classical squared terms of energy contributing to reaction 
absolute temperature 
velocity in x direction 
a species 

distance behind shock 
reactant stoichiometric coefficient 
product stoichiometric coefficient 
mass density of mixture 

standard or sea-level atmospheric density (1.225 x 10 -3 g/cm 3 ) 
constant in equation (10) 

equilibrium value 
forward 
species i 

typical collision partner 
rth reaction 

conditions immediately behind shock 
reference condition at sea level (Earth) 
conditions ahead of shock 


INTRODUCTION 


Round-trip manned lunar and Mars missions may be the next major space exploration goals. For 
these missions, costly propulsive maneuvers will have to be replaced by more economical 
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aeromaneuvers in the upper atmosphere of Earth and the atmosphere of Mars. These aeromaneuvers 
will be performed as high as possible in the atmosphere to minimize heating of the vehicle while 
providing sufficient ambient atmospheric density to generate aeromaneuvering forces. 

Hypervelocity flight in rarefied atmospheres raises questions concerning the thermochemical 
state of the flow field about flight vehicles. Although pioneering study over the past decade has been 
devoted to aeromaneuvering in the upper atmosphere of Earth (ref. 1) and to defining the thermo- 
chemical relaxation phenomena associated with flight speeds related to Earth orbital speeds (ref. 2), 
we must now begin to estimate the extent to which thermochemical relaxation phenomena pertain to 
prospective new missions. These missions (lunar and Mars return) are more energetic than Earth 
orbital missions and involve an atmospheric composition other than the Earth’s (the Mars atmo- 
sphere). Also important is that the vehicles may be larger than those used to date. The estimates 
presented in the current report are especially timely because we are on the threshold of these new 
ventures. 

To estimate the effects of phenomena that are largely beyond our experience, we extend previous 
experimental and analytical efforts, simplifying them to identify gross features. In the current report, 
thermochemical relaxation lengths behind normal shock waves will dominate the considerations. 


AEROMANEUVERING IN EARTH ATMOSPHERE 


Aeromaneuvering at hypervelocity in the upper atmosphere of Earth has been consideied for 
more than two decades, and thermochemical relaxation phenomena in this regime have been studied 
intensely for the past decade. However, the foundation for these studies was laid about three decades 
ago with shock-tube experiments performed by Allen, Rose, and Camm (ref. 3). In those experi- 
ments, a high-speed shock wave was observed as it passed radiometric instrumentation at a station in 
the shock tube. In laboratory coordinates, the radiative intensity behind the moving shock wave 
increased to a peak value and then decayed to an “equilibrium value at some distance behind the 
shock, for each test condition. Spectral measurements were also performed that were useful subse- 
quently for devising molecular models of thermal nonequilibrium. The radiative pulse that peaked 
above an equilibrium value in the experiments was termed “nonequilibrium radiation.” At the time 
the results were only of academic interest, primarily because the dominant heating of the Apollo cap- 
sule occurred at low altitudes where the atmospheric density was higher than the experimental den- 
sity and the flow field about the vehicle was in equilibrium. However, the nonequilibrium radiative 
pulse behind the shock wave in the experiments contained the effects of both thermodynamic and 
chemical nonequilibrium phenomena and represented the coupled thermochemical relaxation behind 
the shock wave— the subject of intense study in recent years. For the current estimates of thermo- 
chemical relaxation effects, the data can be used directly. 

Reference 3, figure 10, plots the product of pressure and the laboratory time after the shock wave 
passes to reach within 10% of thermochemical equilibrium (the ordinate) as a function of shock- 
wave velocity (the abscissa) for a large range of flight conditions in the Earth’s atmosphere. The plot 
can be transformed into the corresponding relaxation length behind a normal shock wave. The result 
appears in figure 1 as three lines near the top center of the figure, labeled 1 cm, 10 cm, and 100 cm. 
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These lines are the loci of the thermochemical relaxation lengths. Some flight trajectories, and other 
information, appear on the figure, to relate the relaxation regime to flight experience. Of particular 
interest is the Aeroassist Flight Experiment (AFE) trajectory. The AFE will examine thermochemical 
relaxation phenomena and their effects, for developing Aeroassist Space Transportation Vehicles 
(ASTV). Such vehicles, which will move space payloads from one orbit to another in Earth-Moon 
space, requiring changes in both orbital altitude and orbital plane inclination, will function in the 
AFE domain. The figure shows that for the AFE at peak heating, the thermochemical relaxation 
length behind a normal shock wave is about 10 cm. That is a significant part of the shock layer 
because the bow shock is about 17 cm in front of the vehicle. Thus a large fraction of the AFE fore- 
body flow field is out of equilibrium thermochemically, which is the intent of the experiment. 

In figure 1 , the point of peak heating on a lunar return vehicle lies below the perigee of the AFE. 
The lunar return vehicle would enter and maneuver in the Earth’s upper atmosphere, similarly to the 
ASTV, and exit the atmosphere to rendezvous with a vehicle in Earth orbit. The figure shows that for 
this vehicle, the thermochemical relaxation length behind a normal shock wave is estimated to be 
between 1 and 10 cm — say about 5 cm. For a large, drag-configured vehicle, the thermochemical 
relaxation length would be of some significance if the shock layer were much thicker than 5 cm, but 
its effects would not dominate the forebody flow field. Thus the importance of the thermochemical 
relaxation phenomena can be estimated from the geometry and size of the prospective lunar return 
vehicle. 

Finally, consider the entry into the Earth’s atmosphere of a manned Mars return vehicle. It may 
perform an aerocapture maneuver as shown in figure 1 , or it may descend to Earth. The entry veloc- 
ity is large, about 14 km/sec. (That speed also includes an emergency abort entry.) The vehicle is 
very energetic, and it experiences peak heating lower in the atmosphere than the other vehicles dis- 
cussed. It is estimated that the thermochemical relaxation distance behind a normal shock wave is 
less than 1 cm, and its effects may be relatively unimportant, depending on the thickness of the for- 
ward shock layer. 


AEROMANEUVERING IN THE MARS ATMOSPHERE 


The composition of the Mars atmosphere is considered to be about 95% CO 2 and 5% N 2 (ref. 4). 
The density of the flight domain is low and thermochemical relaxation should be examined. Because 
the dominant atmospheric species is considered to be CO 2 , it is useful to note some features of the 
basic global conservation equations for compressible flow of CO 2 along a streamline where transport 
phenomena are neglected, but relaxation phenomena are allowed. The conservation equations for 
mass, momentum, and energy for the steady state are, respectively. 


du dp 
p — + u— = 0 
K dx dx 


( 1 ) 


du dp 
pu — = — - 
dx dx 


( 2 ) 
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du dh 
u — + — = 
dx dx 


( 3 ) 


In the absence of relaxation processes, a solution of these equations is that all variables velocity 
(u), density (p), pressure (p), and static enthalpy (h)— are constant. However, for a relaxing flow 
field behind a strong normal shock wave in which the density behind the shock is much larger than 
that ahead of the shock, it can be shown from these equations which variables are most sensitive to 
any process that occurs in the flow; density and velocity can be affected significantly, but pressure 
and static enthalpy are affected very little. This can be shown without solving the equations or speci- 
fying the process. The process enters the problem through the equation of state and whatever sup- 
plemental equations describe the process. Also, since pressure is relatively invariant and density is 
potentially highly variable, temperature is highly variable because of the equation of state. Although 
this argument is simple, it is not obvious or intuitive, and examples that illustrate the sensitivity will 
be provided in the following sections. The concept is useful because before solving the equations, 
the relative behavior of the global flow and thermodynamic variables can be anticipated. This knowl- 
edge may affect the formulation of the method of solving the equations, the devising of numerical 
algorithms, and perhaps the determination of whether a solution can be obtained at all. Thus, for the 
current approach, we know that it is not necessary to deal with large excursions in pressure or static 
enthalpy caused by relaxation processes. 


Thermal Relaxation 

In an early analysis by Howe and Sheaffer of thermochemically relaxing flow in a predominantly 
CO 2 atmosphere, each species was allowed to assume various levels of thermal excitation parametri- 
cally, with no mechanism for thermal relaxation (ref. 5). The analysis showed that the effects of 
thermal excitation on both the flow field and the distribution of chemical species was insignificant. 
(At the time, the study in ref. 5 was not specifically related to Mars entry because the Mars atmo- 
sphere was estimated to contain less than 10% CO 2 , as noted in refs. 6 and 7). 

A subsequent analysis by Hindelang of thermochemical relaxation behind strong shock waves in 
CO 2 showed that the effects of vibrational relaxation and vibration-dissociation coupling are much 
stronger in CO 2 , which has three different vibrational modes, than in diatomic gases, which have 
only a single mode (ref. 8). However, at shock speeds greater than 6 km/sec, the process of finite- 
rate dissociation of CO becomes dominant, and the region of CO 2 relaxation is only a small part of 
the overall relaxation region. Thermally, vibrational relaxation does not have a significant effect on 
the flow field at these higher speeds. 

Recently, Candler obtained a detailed two-dimensional solution for the forebody flow fields of 
vehicles flying at speeds of 6.1 to 7.4 km/sec in the Mars atmosphere, in ambient densities of the 
order of 1(H kg/m 3 (ref. 9). Thermally, two temperatures (translational and vibrational) were con- 
sidered; chemically, 12 reactions and their “reverse” reactions were considered. Thermochemically, 
the flow field was allowed to relax at a finite rate. Candler concluded that “for the cases studied, 
there is very little thermal nonequilibrium due to the very fast vibrational relaxation of CO 2 . 
Although this analysis considered the elements C, O, and N, it confirmed the results of the early 
analyses which were limited to C and O (refs. 5 and 8). 
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For estimating relaxation lengths, the balance of this study will be limited to chemical relaxation 
in a predominantly CO 2 atmosphere. 


Chemical Relaxation 

To proceed with estimates of the chemical relaxation lengths behind strong normal shock waves 
in CO 2 , we write the equation of state for the mixture of CO 2 and its components as 

k 

p = pRT^ n j (4) 

i=l 

where the summation is over all involved species. It is assumed, because of the analyses described in 
the preceding section, that chemical relaxation is characterized by the translational temperature. The 
static enthalpy of each species is 


T 

h i= hP+ f c dT (5) 

JO 1 

and the static enthalpy of the relaxing mixture is 

k 

h(p,p,nj,...,n k ) = ^ njhj (6) 

i=l 

For each reaction r that involves species i, the chemical reactions that pertain to the relaxing flow 
field can be expressed as 

k k k 

X «r,i x i — ^ X P '-i X i (7) 

i=l i=l 

where the arrow indicates that the reaction proceeds to the right at a rate characterized by kf r (It is 
possible, but not generally true, that the reverse reaction is also allowed; for example, the oxidation 
of carbon monoxide in a single step is unlikely in this context.) 

For each species, the variation along the streamline can be written by manipulating equa- 
tions (l)-(7) to form 


d n i _ 1 
dx pu 


I 


(P r,i ^r,i^f 


flipi,)"' 


( 8 ) 
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Equation (8) shows that the relaxation process is coupled to the flow field, and the analysis should 
include flow-field variations that interact with the relaxation process. 

For current purposes, we can allow c Pj /R to be constant for each species. The following expres- 
sion can then be derived from equations (1>— (3) and (5)-(6) (see ref. 10): 
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(9) 


(This differential equation was modified in reference 8 to account for vibrational relaxation.) Equa- 
tion (9) is solved for chemical relaxation behind a strong normal shock wave by use of equation (8), 
which is solved for each species that is participating in the rth reaction at a rate given by 

k fr = 2[PAd 2 /a(s-l)!](2KRT/M*) 1/2 (E fr /RT) s " 1 e" (Efr/RT) (10) 

This form of the reaction coefficient highlights several physical quantities that need to be determined 
(or estimated): P, the sterric factor; s, which is related to the number of contributing energy modes; 
and Ef r , the activation energy for the reaction. For a gas that is relaxing thermally and is character- 
ized by more than one temperature, this expression is difficult to evaluate meaningfully. For a single 
temperature, the behavior of the expression can be examined parametrically, and judgments can be 
made that have some physical meaning. This was done for CO 2 and its components at a time when 
no experimental results were available to the author (refs. 5 and 10), and the results are compared in 
figure 2 with experimental results. For current purposes, only density ratios of the order of p°o/po ~ 
10- 4 are of interest. The flight speeds of interest are from 6 to 9 km/sec; these two speeds are shown 
at the top of the figure with their corresponding temperature ranges, to orient the reader. Two 
reactions are shown in the figure: 


CO 2 + M -> CO + O + M (11) 

CO + M -» C + O + M (12) 

where M is any typical collision partner. The second reaction, equation (12), may be a global reac- 
tion. Its steps were reported by Fairbaim to include the formation and dissociation of C 2 (ref. 1 1), 
although for the elevated temperatures under consideration, that is a rather surprising concept. How- 
ever, the rate coefficients used in the current report are shown by the solid lines in the figure. The 
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rate coefficient for reaction (11) is compared in the figure with the experimental data of Brabbs, 
Belles, and Zlatarich (ref. 12) at a temperature below the range of interest. The predicted rate coeffi- 
cients for both reactions are also compared with the shock-tube data of Davies (refs. 13 and 14). At 
the temperatures of interest, the estimate of the CO dissociation rate is close to the measurements of 
Davies. The estimate of the CO 2 dissociation rate is higher than Davies’ result at about 1 1,000 K and 
is close to but lower than Davies’ result at about 6000 K. In the following analysis, the estimates of 
the forward-reaction rates will be used to estimate relaxation lengths. Generally, the relaxation 
lengths will be close to but less than those that would be predicted by use of Davies’ data. 

A solution for the chemically reacting flow field in CO 2 is shown in figure 3 for a flight speed of 
6 km/sec at ambient density ratio poo/po = 10 -4 (adapted from ref. 10). The distance behind the 
shock is measured along the abscissa, with x = 0 representing the shock location. The distance scale 
is broken in order to capture the variations just behind the shock wave as well as those far down- 
stream of the shock. (A log scale was not used because of the singularity at x = 0.) On the right side 
of figure 3(a), the ratios of equilibrium values to values just behind the shock are shown for the vari- 
ables p, p, u, and T. Density, velocity, and temperature vary significantly in the flow field, as pre- 
dicted previously in the current paper, and pressure and static enthalpy do not. Figure 3(b) shows that 
for this flight condition, only the dissociation of CO 2 to CO and O is important. The dashed curves 
represent an analytic solution wherein the chemistry is uncoupled from the flow field, as described in 
reference 10. 

Figure 4 displays the solution for a higher flight speed: 9 km/sec. The dissociation of both CO 2 
and CO are important at this speed, which is more than twice as energetic as the previous example, 
and the analytic approximation does not represent the chemical variation as well as at 6 km/sec. 

At 9 km/sec, the flow field is sufficiently energetic to partially ionize the reacting mixture; a 
solution that includes ionization is shown in figure 5 (which was adapted from ref. 5). In that solu- 
tion, five reactions were included, some of which were ionizing reactions. The ionizing-reaction-rate 
coefficients were also estimated, but are not assessed for their validity in the current paper. The 
variations caused by ionization can be seen by comparing figure 5 with figure 4. It may be noted in 
figure 5 that ionized species sometimes overshoot their equilibrium values in the hot part of the 
relaxing flow field, as do other species. These overshoots can affect the gaseous radiation and other 
phenomena. 

Figure 6 shows the combined results of many solutions to the problem of chemical relaxation, in 
terms of relaxation length. In this figure the percentage of the total molar change required to achieve 
chemical equilibrium is plotted as a function of distance (or length) behind the normal shock wave. 
Three curves are shown, one for each of three flight velocities, at ambient density pcw/po = 10 -4 . The 
curve for 6 km/sec contains primarily the effects of CO 2 dissociation (see fig. 3(b)). We may infer 
from figure 5 that the curve for 8 km/sec reflects the dissociation of both CO 2 and CO, and a very 
small degree of ionization (ionization-solution results from ref. 5 were included). The 10 km/sec 
curve also reflects the dissociation of both CO 2 and CO, and the degree of ionization for this curve 
was significant. However, the effect of ionization on chemical relaxation length was insignificant, 
possibly because the reactions proceeded more rapidly at this more energetic flight speed (within the 
limitations of this approximate formulation). Note that the curve for 6 km/sec crosses that for 
8 km/sec. For progressively lower flight speeds, the crossing may be expected to progress leftward, 
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and in the limiting case of compressible flow that is not reacting chemically, the curves would cross 
at a relaxation length of zero. We must not speculate too freely, however, because the problem is 
coupled and nonlinear, and the inputs are uncertain estimates. Suffice it to say that the curves cross. 
Also note that for the flight conditions and for all flight speeds shown in figure 6, 90% of the chemi- 
cal relaxation is complete on the order of 10 cm behind the shock wave. 

A careful examination of results in references 5 and 10 leads to the plot in figure 7, which shows 
the chemical relaxation length as a function of flight speed at Poo/po = 10 -4 . Here, relaxation length 
corresponds to the location where 90% of the chemical relaxation is achieved. The relaxation length 
peaks at about 54 cm at a flight speed of about 7.5 km/sec. M. E. Tauber calculated peak heating at 
7.7 km/sec at density ratio p^po = 1.24 x 1(H (private communication, Ames Research Center, 
1990). To the extent that binary scaling applies, the chemical relaxation length scaled from figure 7 
would be 54 cm/1.24 = 43 cm. 

The estimates of the chemical relaxation length behind a normal shock wave in CO 2 at an ambi- 
ent density of the order of Poo/po = 10~^ at flight speeds of about 7 km/sec are dominated by the dis- 
sociation of CO 2 and CO. The estimates are probably meaningful because dissociation rates are 
available. The rate used in the analysis for CO dissociation is close to the measured rate, while that 
for CO 2 dissociation is near or higher than the measured rate. Although the latter may make the 
estimates conservative, the CO 2 dissociates so rapidly that the estimate of the overall length is prob- 
ably affected very little. In any event, these estimates can be refined when more accurate information 
is acquired or when other refinements or judgments are made. 

Finally, shock-tube measurements related to nonequilibrium flow behind normal shock waves in 
CO 2 for a hypothetical high-speed entry into the Venus atmosphere were made by Nealy (ref. 15). 
The nonequilibrium radiative-intensity waveform behind the shock wave was recorded at a wave- 
length of 158 nm. The time for the intensity to peak and then decay to half the difference between 
peak and equilibrium levels was measured. C. Park scaled a corresponding partial relaxation length 
from this data; for a shock velocity of 10 km/sec at Poo/po = 5 x 10~^, it was about 12.5 cm (private 
communication, Ames Research Center, 1990). Binary scaling for poo/po = 10 -4 yields a partial 
relaxation length of 5 x 12.5 cm = 62.5 cm. Although a rigorous comparison with figures 6 and 7 is 
not feasible, the relaxation length inferred from the experiments is large — probably larger than the 
present estimate at that high-speed condition, and perhaps two or three times that value. The conclu- 
sion that nonequilibrium flow is significant and should be considered in Mars flow-field analyses is 
thus reinforced by both analysis and experiment. 


CONCLUDING REMARKS 


Estimates of thermochemical relaxation lengths behind normal shock waves in air and CO 2 were 
made, to assess the significance of nonequilibrium phenomena in studies of forebody flow fields for 
hypervelocity aeromaneuvering in the atmospheres of Earth and Mars. For Earth the thermochemical 
relaxation length is a large fraction of the bow-shock standoff distance for the Aeroassist Flight 
Experiment currently being planned. That is desirable because one of the objectives of the experi- 
ment is to study the nonequilibrium phenomena. For manned return from the Moon to Earth by use 
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of aeromaneuvering vehicles, thermochemical relaxation may be insignificant near peak heating if 
the vehicle configuration is large and blunt. For the manned Mars return aerocapture maneuver, 
thermochemical relaxation near peak heating seems to be of little importance. 

Thermodynamic relaxation for entry into the Mars atmosphere seems to be of little importance 
on the vehicle forebody for flight conditions of current interest. On the other hand, the chemical 
relaxation length seems to be significant. It is estimated to be about half a meter at peak heating. This 
result can be scaled to other flight conditions. 
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Figure. 1 The Earth aerothermodynamic flight regime, with the thermochemical relaxation domain 
highlighted. 
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Figure 2. Forward-reaction-rate coefficients. 
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Figure 6. Relaxation length behind normal shock wave in CO 2 — flow field coupled to chemistry; 
poo/po = 10- 4 . 



Figure 7. Chemical relaxation length (for 90% relaxation) behind normal shock wave in CO 2 ; 

poo/po = 10" 4 . 
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